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1 Introduction 

The Method of Moments (MoM) introduced in the previous lecture is widely used for solving 
integral equations in the frequency domain. Yet, some attempts have been made recently at the 
use of the MoM in the time domain. We shall briefly expose this approach here. 

2 Time domain equations 

The first step is of course to write Maxwell's equation and all other relations (constitutive 
relations and continuity) in time domain: 

VxE(f,t) = -^B(f,t)-M(f,t), V x H(f,t) = ^D(f,t) + J(f,t) , (la) 

V ■ B(f,t) = m(f,t) , V ■ D(f,t) = p(f,t) , (lb) 
D(f,t) = eE(f,t), B(f,t) = fiH(f,t), (lc) 

V ■ J(f,t) + —p(f,t) = 0, V • M(f,t) + — m(f,t) = 0. (Id) 

For the time-domain MoM, it is easier to work with the potentials, and make use of the 
well-known "retarded potentials" theory. In view of doing this, we write the definition: 

H(f, t) = — V x A(f, t) (2a) 

E(f,t) = -V<t>(f,t)--A(f,t). (2b) 

Both the vector potential A and the scalar potential (f> satisfy the wave equation which, in 
time-domain domain, writes: 

d 2 

V 2 A{f, t) - eoHo-T^Mr, t) = -H>J(r, t) , (3a) 
V 2 0(f,t) - e ^<t>(f,t) = . (3b) 

These potentials are linked by the time-domain Lorentz gauge: 

V-A(f,t) + eono^<Kr,t) = 0. (4) 
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Section 2. Time domain equations 



We can denned also a time-domain Green's function which satisfies the time-domain scalar 
equation: 



( V2 " ^^9{r , f, t, f) = -S(r - f) 6(t - t') , 



which solution is (in free-space): 

g(f,f',t,t') = { 



i 6 (t-i/-\^fl) t>t>, 



4-7r|r— r'| 





t < t'. 



From this, the solution to the wave equation for A and <p can be written as: 

, J(r', t - R/c) 



A{f,t)=no / dv' / dt'J(f',t')g(f,f',t,t') = no / dv' 

JV J~oo JV 



AttR 



(5) 



(6) 



(7a) 



<p(f, t) = dv 
JV 



, p(r,t- R/c) 



4tt€ R 



(7b) 



where R = \f — f'\. These wave equations are known as the time retarded potentials, and 
essentially say that the potential (either A or 4>) can be calculated at a given point in space f 
and given time t from all previous times. 

From these equations, we can calculate the space-time electromagnetic fields: 



H(f,t) = j- f d^'Vx^, r = t-R/c, 



E(f,t) = --±- / d,'V%)-f Id, 



47reo J v R 
Let us continue with the electric field first: 



fV7 p{f,T) [ dv , d J{?,t) 

4tt ./ dt R 



E(f,t) 



1 



dv' 



^V P (f',T) + p(f',r)V- 



47re j v 

At this point, we need to use the following relations: 

R 



Mo 

47T 



VR 
1 



V 

Vp(f',T 



R' 
_ R 
R ~ ~W ' 



p (f', T Vr = - Vi?— p f',r = -- -— p f,r 
ar c or c Rot 



We can therefore continue with the electric field as: 



1 R d . _. R ._. 

C^cV P(f ' T) + i?3 P(f ' r 



/'° / ,1 ' 1 ,V 7/-/ \ 



4tt€o 



dv' 



^ p(f '' T) + ^ (f, ' T) 



4"7T 

Ri 4tt / Rdr [r,T) - 



(8a) 
(8b) 

(9) 

(10a) 
(10b) 
(10c) 



(11) 
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We can perform the same type of calculations for the magnetic field using the relation 



We get: 



H(r, t) = i- I dv' 



1 R d - f R - f 

C^ X ^ J(r ' T) -i?3 Xj(r ' r) 



(12) 



(13) 



Upon gathering the expressions for the electric and magnetic field, we eventually get: 

R ixq d 



E(f, t) = — [ dv' 
4vr J 

H(f, t) = -J- [ dv' 
4vr J 



^J(f'.r) + Ij(f'.r) 



e R 2 R dr 
R 



J(f',r) 



(14a) 
(14b) 



Upon using the boundary conditions for the electric and magnetic field, we construct the 
integral equations in a standard way: 



• EFIE: hx{E i + E scat ) = on PEC surface 

%t) + ^x J ds' [...] 



n x 



(15) 



• MFIE: h x (H i + H scat ) = J s . 

As we have seen before (in a previous class), this integral equation is expressed in terms 
of the principal value of the integral with a 1/2 additional factor. Thus: 



- J(f, t) = hx lP(f, t) + -*-n x PV [ ds' [. . .] 
2 47T J 



(16) 



For the sake of comparison, we can write the MFIE in the frequency domain and in the time 
domain: 



J(f) = 2n x H\f) + 2h x PV j ds' J(f') x V'g(f, f') 



feS. 



J(f,t) = 2nx H\f) + — h x PV / ds' 

27T 



^J(f'.r) + ij(f'.r) 



i?2 



(17a) 
(17b) 



Note that in the principal value, we essentially exclude the part for which R = 0. Since 
r = t — R/c and R ^ 0, we always have that t < t. The time domain equations therefore 
state that the current at location f and time t is equal to a known term 2h x H l (f,t) plus a 
term (integral) known from the past history of J. This is the basis for solving the time domain 
integral equation by iterative methods, the most well-known one being the marching- on-in-time. 
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Section 3. The marching-on- in-time technique 



3 The marching-on-in-time technique 

3.1 General equations 

The integral equation can often be cast in the following form: 

J(f,t) = J i (f,t)+ I dv' [ dt'T(f,r,t-t')-J(r,t'). (18) ^qTTq 

Js Jo 

l eg ,10 - 

Note that we have a time integral also as in Eq. (ttojj J(r , t ) has not yet been set to satisfy 
any causality condition. Hence, we must then impose J l (f,t) = for t < 0, f G S. 

In order to apply the MoM, we discretize the current both in space and in time: 

M N 
m'=l n'=0 

where P denotes the simple pulse function. In addition, we also apply point-matching, which 
means that we take the following testing functions: 

W mn (f m , t n ) = S(f - f m ) S(t - nAt) = S(f - f m ) 5(t - t n ) , (20) 

where we take At = mm{R mm //c}, R m m> = \r m — f' m \. The method is best illustrated on the 
following example. 

3.2 Example 

Let us consider a ID example governed by the following integral equation: 

/ x ° \x — x'\ 
K(x,x') f(x',T)dx' , xe[-x ,x ], t = r(x,x',t) = t . (21) 
-x C 

Let us chose the following expansion for /: 

N J 

f(x', t) ~ 2 j'tf, t) , (22) 

i'=i j'=i 



where the pulse basis functions are defined as 



Pi'f = < 



1 for x' e [x v -%,Xi> + 4*-] and t € [tj> - f,t f + f ] 
elsewhere. 



(23) 



Note that we use the definitions: X{> = i'd x , tj> = j'dt, and d x = cdt- In order to apply point 
matching, we take the following testing functions: 

Wij(x, t) = S(x - Xi) S(t - tj) . (24) 
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Upon expanding and testing, we get: 



N J 



/xo 
K(xi,x') ^2 /2 o.i'j'Pi'j'(x',T)6(t - tj) 
-x v=x , /=1 



i'=l 



/ Ax ai'j/Pi'j'(x' ,T)K(id x ,x')S(t — tj) 



Coming back to the definition of r, we write (with the test and the expansion): 



t = tj = jd t - \i - i \ — = (J - \i - i \)d t , 



such that the coefficient a^ji becomes The integral equation becomes: 



N 



9ij ~ ^2 a i',j-\i-i'\ 
i'=l 



(i'+h)d x 



(i'-\)d x 



dx'K(id x , x') . 



We can define the term 



r(i'+h)d* 
I dx'K(id x ,x') 

J{i'-\)dx 



and rewrite the previous system as 

N 

9ij = ^ Z ii' ai\j-\i-i'\ 



i'=l 



Zii CLij + Zij-i di-ij-i + Zij-2 Oi-2,j-2 + • • • + Zn a± t j-i + i 

+ Zi : i + i Oi+ij-i + Zi y i + 2 Qi + 2,j-2 + • • • + Zi^n a>N,j-N+i ■ 



(25) 



(26) 



(27) 



(28) 



(29) 



In this equation, only the first term involves time step j, all the others terms being at j — 1, 
j — 2, ... Therefore, we can solve for a^: 



1 



Hj 



N 



9ij — ^2 a i',j-\i-i'\ 



i'=l 



(30) 



The value of all are known for k < j, so that is completely specified in closed form by 
those and the present value of gij. This process is known as a ID march-on-in-time approach. 

Time-domain MoM is nowadays in its early stage and, although it has been successfully 
applied to various simple situations, still suffers from numerical instabilities. More work is in 
progress... 



